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Abstract 

We study the polarization transport of transverse phonons by adopting a new approach based 
on the quantum mechanics of spin-orbit interactions. This approach has the advantage of being 
apt for incorporating fluctuations in the system. The formalism gives rise to Berry effect terms 
manifested as the Rytov polarization rotation law and the polarization-dependent Hall effect. We 
derive the distribution of the Rytov rotation angle in the presence of thermal noise and show that 
the rotation angle is robust against fluctuations. 
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I. INTRODUCTION 



Berry phase is a non-integrable phase factor arising from the adiabatic transport of a 



system around a closed path in its parameter space 



ij. Geometrically, it originates from 



parallel transport in the presence of a gauge connection in the parameter space m. This 
phase factor is amazingly universal and appears in many different contexts sl- Here, we 
explore some of its physical manifestations in acoustics. 

The analogy between the linearized equations of elasticity and Maxwell equations is well- 
known (see e.g. ^). One can take advantage of this analogy and predict new phenomena 
for acoustic waves from their optical counterparts. In particular, polarization phenomena in 
optics can be mapped onto transverse acoustic waves. For example, the Rytov polarization 



rotation law [5] and the polarization-dependent Hal 




effect in optics, which are manifestations 



I— — 1 

of the Berry effect for photons [6|, [7| 

acoustic waves as well Q, Q, ^ 

Acoustical Berry effects in inhomogeneous media were first established 13] via the post- 

S metric acoustics approximation (the acoustical analogue of the post-geometric optics 
, Q). We have recently treated such media, and in particular those with periodic den- 
sity variations (phononic crystals), via another semiclassical approach (the wave packet 
approximation) [20]. A complication that arises in acoustical polarization transport in in- 
homogeneous media is the coupling between the transverse and longitudinal waves caused 
by the spatial variations of density and/or elastic coefficients. This means that polarization 
cannot be maintained. However, the coupling vanishes in the adiabatic limit of small spatial 



variations which renders polarization an adiabatic invariant [19|, |20]. In the present study, 
we introduce yet another approach that is based on the quantum mechanics of spin-orbit 
interactions. This approach has the advantage of being apt for incorporating fluctuations 
in the system. The formalism gives rise to a Berry phase describing the rotation of the 
polarization vector (the Rytov law), and a Berry curvature in the semiclassical equations of 
motion deflecting the phonons depending on their polarization (the polarization dependent 
Hall effect). These, of course, reproduce the results of [l^ (and 20], as far as non-periodic 
inhomogeneous media are concerned). We analyze the effect of classical noise on the Rytov 
law and derive the distribution of the Rytov rotation angle in the presence of thermal noise. 
It is found that the latter is robust against ubiquitous thermal fluctuations in the system. 
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II. THE QUANTUM MECHANICAL APPROACH 



The dynamics of the displacement vector filed U(x, t) in an elastic medium is governed 
by (see e.g. 0) 



(summation convention implied) where p is the mass density of the medium and cTjj is the 
stress tensor. In an isotropic medium the latter takes the form 

aij = X6ij dkUk + p{diUj + djUi) 

A and p being the Lame coefficients. We consider elastic waves in an isotropic medium whose 
density and/or Lame coefficients vary with position. By decomposing U into a solenoidal 
(transverse) and an irrotational (longitudinal) part, one can show 2l|] that the shear trans- 
verse and the compression longitudinal modes propagate independently in a strictly homo- 
geneous medium. Introducing sufficiently small (adiabatic) spatial variations, the modes 
can still be regarded as independent to arbitrary accuracy [lol, 3 • The adiabatic approxi- 
mation eliminates the coupling between the transverse and longitudinal modes and renders 
polarization an adiabatic invariant. 

The transverse mode constitutes a two- (polarization) state system. Each state is involved 
in a double degeneracy in the absence of inhomogeneity. The double degeneracy is the po- 
larization degeneracy; transverse waves with different polarizations have the same dispersion 
in a homogeneous isotropic medium 22|. In inhomogeneous media, the refractive- index gra- 



dient lifts this polarization degeneracy by coupling the polarization and the translational 



degrees of freedom 



19| . The coupling also affects the direction of propagation of the phonon, 



which is given by its momentum p. Such spin-orbit interactions also play significant role in 



the spin transport of photons and electrons [12|, ll3|, l23|, l2j, l25]. The ensuing Berry effects 
are simply manifestations of this interaction. This motivates an alternative approach to 
the problem that relies only on the quantum mechanics of the spin-orbit interaction. Thus, 
following Berry [Ij, we can derive the geometric phase: The parameter space for adiabatic 
excursions, here, is the momentum space. The spin-orbit Hamiltonian can be represented 
by a (2 X 2) Hermitian matrix coupling the two polarization states. Let us take the point 
in the parameter space at which the sates are degenerate, as the origin. With reference to 
this point, the Hamiltonian can be expanded to first order in p. The most general such 
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matrix satisfying the given conditions depends on the components of p, and by hnear trans- 
formation in the momentum space can be brought to the form CTjpj, where (jj are the Pauh 
matrices. The two hehcity eigenstates |cr, p), with eigenvalues ap (cr = ±1,]? = |p|), form 



the non-degenerate normal modes, exactly as for photons 
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231]. These intersect conically 



at the origin (the degeneracy point). The Berry curvature is calculated [Jj to be 

Vp X ((T,p| iVp|o-,p) = -ap^^p 

which is the field of a magnetic monopole of charge —a situated at the origin of the momen- 
tum space. It corresponds (in appropriate gauge) to the Berry connection 

(cr,p| zVp|o-,p) = crA_L 

where A^(p) = j^q;t^^(~P25Pi5 0). The resulting geometric Berry phase a/c-rfp ■ A_l, 
where C is the phonon trajectory in momentum space, is of opposite signs for the two 
polarizations. Therefore, for a linearly polarized wave, this Berry phase leads to the rotation 
of the polarization plane through the angle 

7 = / ■ (ip = / cos 9 dip (1) 

Jc Jc 

0{(f) being the zenith (azimuth) angle in the spherical polar coordinates of the momentum 
space. This is the Rytov law for transverse phonons. 

Because the spin-orbit interaction introduces a gauge potential (Berry connection) in the 
momentum space, the position operator of a transverse phonon of helicity ah acquires an 



anomalous contribution according to x — > r = x — ahA± 2^. This is in perfect analogy 
with the electromagnetic interaction, with the role of the position and momentum inter- 
changed and helicity replaced by charge, of course. (In the post geometric optics/acoustics 
approximation, the anomalous contribution stems from the representation in which the 
wave equation is diagonal 12, Q, [l^.) The physical (observable) position coordinates are 
no longer the canonical coordinates x = ihVp, but are now r = hliVp — aA±), which are 
non-commutative: 

r 1 • i;2 Pk 

[ri,rj\ = lah Sijk^. 

While \pi,Pj] = and [^'i,^^] = ih6ij, this provides an example of non-commutative quantum 
mechanics. Using these commutation relations, the Heisenberg equations 



yield the semiclassical equations of motion 

p = -VrH, r = VpH + ah 



P X P 

p3 



(2) 



to first order in h. These, of course, reduce to the standard ray equations of geometric 
acoustics in the classical limit h —>■ 0. In this limit, helicity vanishes and the right/left 
circularly polarized waves follow the same trajectory. However, according to ([2]), these rays 
now split due to the effect of the Berry curvature of the momentum space. The deflections 
from their classical (geometric acoustic) trajectories are given by 

p X dp 



6r 



a 



h 



(3) 



ic p-" 

The resulting displacements are, therefore, locally orthogonal to the directions of motion of 
the phonons. This, which is a general feature of spin transport 1^, [ 
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30|, 



constitutes the polarization-dependent Hall effect of transverse phonons. In the present 
formulation, the effect is a direct consequence of the non-commutativity of the physical 
coordinates resulting from the spin-orbit interaction. 

The above findings for Berry effects, equations ([I]) and ([3]), coincide with the results of 



IS 



201], which were obtained differently. 



III. RYTOV LAW IN THE PRESENCE OF CLASSICAL NOISE 

Fluctuations in the system may directly affect the refractive index by randomly perturb- 
ing the density and/or the elastic coefficients of the medium. An ubiquitous example is 
provided by thermal fluctuations. As such, they cause the direction of phonon propagation 
to fluctuate, so that 

p(t) =Po(t)+N(t) 

where suffix indicates the absence of noise and N(t) is a noise term- a random process with 
zero average and small amplitude compared to po. Then, the spin-orbit Hamiltonian aiPi 
still leads to the result ([1]), but the trajectory C is now a fluctuating trajectory- it fluctuates 
about the noiseless evolution trajectory (assumed cyclic with period T). Thus, using ([1]), 
the resulting change in the Rytov rotation angle 7 during time T is given by 

A7(T) = ^ /^(cos^ - coseo)dt (4) 
T JO 
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In case, by virtue of the random noise, p does not return to its original direction, a non- 
cyclic contribution also appears. This term must be removed according to the definition of 
Berry phase for non-cyclic evolution 3l|, so that the above result still holds 32|]. Writing 
the integrand in (jl]) as 

cos 6 — cos Uq = 

P Po 

and expanding p in terms of po to first order in the noise, yields the result 

A7(T) = ^ r - p„ . N)*. (5) 

T Jo Po Po 

This is the Rytov law in the presence of noise. Thus, (A7) = 0, which means that the average 
value of the Rytov angle coincides with its noiseless value. To characterize the probability 
distribution of A7 further, we need a definite model for the noise. The uncorrelated noise 
defined by 

(iV,(t)) = 0, {N,{t)N^{t')) = 2D6{t - 

where angular bracket denotes ensemble average, serves suitably in view of the physical 
nature of thermal fiuctuations. Then, yields 



for the variance of the distribution, where bar denotes average over the evolution time. The 
effect of fiuctuations, thus, diminishes as 1/T so that 7 coincides with its noiseless value 
in the adiabatic limit T ^ 00. In other words, the Rytov rotation angle is robust against 
ubiquitous thermal fiuctuations in the system. 
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